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The relationship between the densities of ground-state wave functions (i.e., the minimizers of 
the Rayleigh-Ritz variation principle) and the ground-state densities in density-functional theory 
(i.e., the minimizers of the Hohenberg-Kohn variation principle) is studied within the framework of 
convex conjugation, in a generic setting covering molecular systems, solid-state systems, and more. 
Having introduced admissible density functionals as functionals that produce the exact ground 
ground-state energy for a given external potential by minimizing over densities in the Hohenberg- 
Kohn variation principle, necessary sufficient conditions on such functionals are established to ensure 
that the Rayleigh-Ritz ground-state densities and the Hohenberg-Kohn ground-state densities are 
identical. We apply the results to molecular systems in the Born-Oppenheimer approximation. For 
any given potential v £ L 3/,2 (R 3 ) + L°°(R 3 ), we establish a one-to-one correspondence between the 
mixed ground-state densities of the Rayleigh-Ritz variation principle and the mixed ground-state 
densities of the Hohenberg-Kohn variation principle when the Lieb density-matrix constrained- 
search universal density functional is taken as the admissible functional. A similar one-to-one 
correspondence is established between the pure ground-state densities of the Rayleigh Ritz variation 
principle and the pure ground-state densities obtained using the Hohenberg- Kohn variation principle 
with the Levy-Lieb pure-state constrained-search functional. In other words, all physical ground- 
state densities (pure or mixed) are recovered with these functionals and no false densities (i.e., 
minimizing densities that are not physical) exist. The importance of topology (i.e., choice of Banach 
space of densities and potentials) is emphasized and illustrated. The relevance of these results for 
current-density-functional theory is examined. 
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I. INTRODUCTION 


The concept of modeling an electron gas and more generally all electronic systems using only the density 
p goes back to the early days of quantum mechanics and the independent work of Tomas and Fermi in 1927 
1, 2]. However, it was only with the seminal work of Hohenberg and Kohn in 1964 that density-functional 
theory (DFT) was shown to be—in principle, at least—an exact theory 0. In 1983, Lieb formulated DFT 
rigorously for electronic systems in R 3 , using concepts from convex analysis 0 . Today, DFT is the most 
popular computational method for many-electron systems in chemistry and solid-state physics. 

In an external electrostatic scalar potential v, the ground-state energy E(v) of an iV-electron system is 
obtained from the Rayleigh-Ritz variation principle 

E{v) :=in£(if\H(v)\if), (1) 

where H(v) is the Hamiltonian and the infimum extends over all L 2 -normalized wave functions for which 
the expectation value makes sense. With each minimizing ground-state wave function if in Eq. m, there is 
an associated ground-state density p. In DFT, the ground-state energy E(v) is instead obtained from the 
Hohenberg-Kohn variation principle , 

E(v) = inf(FoO) + (w|p)), (2) 

p 

where Fq ( p ) is the universal functional and the minimization is over all densities for which the interaction 
(v\p) = /u(r)p(r)dr is meaningful. Note that we have not specified the spatial domain nor the nature of the 
external potential v. Thus, the minimization problems mm are so far not fully defined from a mathematical 
perspective. Moreover, the functional Fq is not unique—any functional Fq that gives the correct ground- 
state energy in the Hohenberg-Kohn variation principle is said to be ‘admissible’. Common expositions 
of DFT use the Levy-Lieb constrained-search functional Fll 3,1, the density-matrix constrained-search 
functional Fdm 0 , and the Lieb functional F, all giving the same ground-state energy E(v) for a large 
class of potentials. Specializing to electrons in R 3 and assuming that the external potential v is sufficiently 
well-behaved, these functionals are given by the following expressions: 

*ll(p) := mf{(if\H 0 \if) 1€ ifg(R 3JV ), \\ifh = 1, -> p} , (3) 

^dm(p) := inf|Tr(77 0 r) | T = Efc \k\ifk)(ifkl Efc E = 1 , E > 0, if k £ H^(R 3N ), \\tfkh = 1, if ^ p}, (4) 

F{p) := sup (E(v) - (v\p)). (5) 

vex 1 

Here Hg (R 3;V ) is the first-order Sobolev space Hg (R 3;V ) with proper permutation symmetry due to spin, 
H o = T + W is the sum of the kinetic and inter-electron interaction operators, and if p and T h->- p 
indicate that the wave function if and density matrix T, respectively, have density p. In the definition of 
the Lieb functional, X is a Banach space in which densities are embedded, with X' as its dual, consisting 
of potentials. We note that the Lieb functional is by construction lower semi-continuous and convex, being 
the conjugate function to E in the sense of convex analysis. 

Whereas the density functionals Ell, ^dm and F are all admissible and therefore give the same total 
energy in the Hohenberg-Kohn variation principle, they may have different minimizing densities. In this 
article, we study the relationship between the densities in the Hohenberg-Kohn and Rayleigh- Ritz variation 
principles. In particular, we establish what conditions must be imposed on an admissible density functional 
to ensure that the ground-state densities (i.e., the minimizing densities) in the Hohenberg-Kohn variation 
principle are the same as the densities of the ground-state wave functions (i.e., the minimizing wave functions) 
in the Rayleigh-Ritz variation principle. 

While our treatment covers the standard DFT setting outlined above, it is motivated by a study of 
current-density-functional theory (CDFT), where the iV-electron system is subject to an external magnetic 
vector potential A in addition to the scalar potential v 0-0|, introducing the paramagnetic current density 
j p £ L 1 (R 3 ) as an additional variable in the Hohenberg-Kohn variation principle: 

E(v. A) = inf (F 0 (p, j p ) + (v + ±A 2 | p) + (A | j p )) , (6) 

PiJ P 
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where (A | j p ) = J A(r) j p (r) dr. The CDFT generalization of the Levy-Lieb functional is the Vignale-Rasolt 
functional TVr(p, j P ); there is likewise a generalization of the density-matrix constrained-search functional 
Tdm and the Lieb functional F to include a current dependence. The mathematical properties of these 
functionals are not as well understood as in the standard DFT setting. The results presented in this article 
are in part meant to be a step towards such an understanding. 

The article is structured as follows. In Sec. in we describe DFT from a generic and abstract point of view, 
introducing arbitrary admissible density functionals Fo an d some concepts of convex analysis. We obtain 
several results that characterize the relationship between minimizers of the Rayleigh-Ritz variation principle 
and of the Hohenberg-Kohn variation principle. In Sec. lIIIl we apply our findings to the standard DFT 
of atoms and molecules in the Born-Oppenheimer approximation. The importance of the topology of the 
underlying density space is emphasized. In Sec. lIVI we discuss CDFT, and finally, in Sec. [V] we summarize 
and draw some conclusions. 


II. DFT FROM AN ABSTRACT POINT OF VIEW 
A. Admissible universal density functionals 

Except for the trivial case of one electron, an explicit formula for any of the standard density functionals 
used in the Hohenberg-Kohn variation principle in Eq. m is not known. Moreover, the mathematical 
analysis of Eq. © is difficult without further assumptions. In his work [dj, Lieb placed DFT for electronic 
systems in R 3 on a firm mathematical ground using the language of convex analysis I,® , the natural 
setting for problems such as Eq. ©• The starting point is to embed the densities in a Banach space X and 
to consider potentials in the dual space A', thereby making the interaction (v \ p ) well-defined and continuous 
in both arguments. We then obtain the ground-state energy as a map E : X' —> R U {— 00 } given by 

E(v) = inf (F 0 (p) + (u|p)), V® 6 T', (7) 

p£X 

where Fo : X — > M U {+ 00 } is the universal density functional. Being the pointwise infimum of a collection 
of affine maps of the form v (v | p) + Fo (p ), the ground-state energy in Eq. © is automatically concave 
and upper semi-continuous with respect to the topology on X'. Note that we allow E and Fq to be infinite. 
The effective domain of a function is the set where the function is finite—for example, dom(F?) = {v £ X 1 \ 
E(v) > — 00 }. 

In Ref. [3], Lieb (and Simon) proved several important results in the context where X = F 1 (R 3 ) flF 3 (R 3 ), 
with dual X' = F 3 / 2 (R 3 ) + F°°(R 3 ). However, this choice of X is not unique, and the derived results may 
depend on X. We emphasize that, even if our notation suggests the standard DFT setting, it covers also 
CDFT and other settings. 

A function Fq : X —> RU {+ 00 } is said to be an admissible density functional if, for each potential v £ X ', 
it gives a correct ground-state energy by the Hohenberg-Kohn variation principle in Eq. © — that is, if it 
produces identical results with the Rayleigh Ritz variation principle in Eq. ©. Note that we do not require 
that there exists a minimizing density p for a given potential v, even in those cases where v supports a ground 
state in the Rayleigh Ritz variation principle in Eq. ©. It is sufficient that a minimizing sequence {p n } 
can be found. Clearly, for any pair (p, v) £ X x X\ an admissible density functional and the ground-state 
energy satisfy the Fenchel inequality 


E{v) < F 0 (p) + (v\p). (8) 

We are here interested in characterizing those pairs (p,v) that saturate Fenchel’s inequality, E(y) = Fo(p) + 
(u|p); in other words, those p £ X that are minimizers in Eq. © for a given v £ X'. 

Since the universal functional is in general not differentiable (see Refs. UM for the standard DFT case), 
we cannot write down an Euler equation for the solution of the minimization problem Eq. ©. Moreover, even 
if Fo were differentiable, the Euler equation would in general be a necessary but not sufficient condition for 
a global minimum in Eq. ©. Instead, we use the concept of subdifferentiation to characterize the solution. 
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The admissible density functional Fq is said to be subdifferentiable at p £ X if Fq (jj) £ K and if there exists 
u £ X' , known as a subgradient of Fq at p, such that 

F 0 (p')>F 0 (p) + (u\p'-p), Vp'£X, (9) 

meaning that Fq touches the affine map p’ i->- Fq(p') + (u \ p’ — p) at p and lies nowhere below it. The 
subgradient is thus a generalization of the concept of the slope of a tangent to the graph of Fq at p. The 
subdifferential d~ Fq(p) is the (possibly empty) set of all subgradients of Fq at p: 

d~F 0 (p ) = { u £ X' | Fo(p') > F 0 (p) + 0 u\p ' - p), V £ X, F 0 (p) £ R } . (10) 

This set is convex Note that the subdifferential is the empty set whenever Fo(p) = +oo. Assuming that 
— v is a subgradient of Fq at p , we obtain by simple rearrangements 

Fo(p') > F 0 (p) + (-v | p’ - p) 

Fo(p) + {v | p) < F 0 (p') + (v | p') 

Fo(p) + {v |p) = inf (Fo(p') + {v \p')) = E{v), (11) 

p'GA 

We have thus shown the following sufficient and necessary condition for a minimizing density in the 
Hohenberg-Kohn variation principle in Eq. ©: 

Proposition 1. Let Fq : X ► K U {+oo} be an admissible density functional so that Eq. (Q holds. Let 
v £ X' and p £ X be given. Then, 

E{y) = F 0 (p) + (v\p) -v £ d~F 0 (p). (12) 

B. Lieb’s universal density functional 

From the Fenchel inequality in Eq. ([8]), we obtain by a trivial rearrangement the equivalent inequality 

F 0 (p) >E(v) - (v\p), (13) 

stating that each admissible density functional Fo(p) is an upper bound to p i —> E(v) — (v\p) with respect to 
all variations in v £ X'. We now define the Lieb universal density functional F{p) as the least upper bound 
to E(v) — (v | p) for all v £ X': 


£ d F 0 {p) 


F(p) = sup (E(v) - (v | p)), Vp £ X, (14) 

vex' 

which in the following will be referred as the Lieb variation principle. The Lieb functional is clearly a lower 
bound to all admissible density functionals: 


F(p) < F 0 (p). (15) 

Since the Lieb functional by construction also satisfies the Fenchel inequality in Eq. ©, we obtain: 
Proposition 2. The Lieb functional is an admissible density functional: 

E{y) = inf (F(p) + (v | p)), Vu £ X'. (16) 

pex 

Proof. E(v) < inf pg x (F(p) + (v \ p)) < inf p&x ( F 0 (p ) + (v \ p)) = E(v). □ 

The Lieb functional is related to the ground-state energy in a special, symmetrical manner—compare the 
Lieb variation principle in Eq. (11411 with the Hohenberg-Kohn variation principle Eq. (1161) . In the language of 
convex analysis, E and F are said to be skew conjugate functions jl2j : the function E is concave and upper 
semi-continuous, whereas its skew conjugate function F is convex and lower semi-continuous. Convexity of 
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F and concavity of E mean that, for each pair pi,p 2 G X, each pair v\,v 2 £ X ', and each A G (0,1), we 
have 


F(X Pl + (1 - A ) P2 ) < XF(pi) + (1 - A )F(p 2 ), (17) 

E {Am + (1 - A)u 2 ) > A E( Vl ) + (1 - A )E{v 2 ), (18) 

whereas lower semi-continuity of F and upper semi-continuity of E imply that 

liminf F{p) > F(p 0 ), (19) 

p->Po 

limsuplH(ri) < E(v o). (20) 


These properties follow straightforwardly from Eqs. m and (fl6l) . It is a fundamental result of convex 
analysis that there is a one-to-one correspondence between all lower semi-continuous convex functions on X 
and all upper semi-continuous concave functions on X' SEE EH, the correspondence being as in Eqs. m 
and (flGl) . It follows that the Lieb functional F is not only a lower bound to all admissible density functionals 
but also the only admissible density functional that is lower semi-continuous and convex with respect to the 
topology on X. It is a trivial but important observation that each property and each feature of E are, in 
some manner, exactly reflected in the properties and features of F and vice versa. Note, however, that the 
Lieb functional depends explicitly on X, just like E depends on X'. 

Just like there may happen to be no minimizing density in the Hohenberg-Kohn variation principle, there 
may happen to be no maximizing potential in the Lieb variation principle. To characterize maximizing 
potentials, we introduce superdifferentiability by analogy with subdifferentiability. The ground-state energy 
E is said to be superdifferentiable at v £ X' if there exists an element p G X, known as a supergradient of 
E at v, such that 

E{ v') < E(v) + (v' -v\ p), Vv' G X'. (21) 

The superdifferential d + E{y) is the (possibly empty) convex set of all supergradients of E at v: 

d + E{v) = { p G X | E(v') < E{v) + (p\v' - v), Vv' G X', E(v) G R } . (22) 

In exactly the same manner that we proved Eq. m, we obtain the following necessary and sufficient 
conditions for the existence of a maximizing potential in the Lieb variation principle: 

E(v) = F(p ) + (v\p) <£=> p G d + E{v). (23) 

Note carefully that this result holds only for the Lieb functional, not for an arbitrary admissible density 
functional. Combining Eqs. m and (1231) . we arrive at the following characterization of the optimality 
condition in Eqs. ITU) and HUl) : 

Proposition 3. If F : X M. U {+oo} is the Lieb functional and E : X' —> R U {—oo} is the ground-state 
energy , then 

E(v) = F(p) + (v\p) <*=>■ —v G d~F(p) <£=> p G d + E{v). (24) 

For a general admissible density functional Fq, these conditions are not equivalent: there may exist 
p G d + E{y) such that — v ^ d~F 0 (p) when F 0 F. On the other hand, the converse statement, — v G 
d~F 0 (p) => p G d + E(v), holds for any admissible density functional F 0 . To prove this, assume that 

— v G d~F 0 (p). According to Eq. (fT^l) . we then have F 0 (p) = E{v) — (v \ p). At the same time, the Fenchel 
inequality holds for any admissible density functional Fq: 

E{u) < F 0 (p) + (u |p), VuG X'. (25) 

Substituting Fq(p) = E{v) — (v \ p) in this inequality, we obtain 

E(u) <E(v) + (u-v\p), VuGX', (26) 

implying that p G d + E(v). We have thus proved the following: 
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Proposition 4. If Fq : X —t RU{+oo} is an admissible functional and E : X' —> RU{—00} the ground-state 
energy, then 


E(v) = F 0 (p) +(v\p) <*=>■ -v£d F 0 (p) => p £ d + E(v). ( 27 ) 

In general, we have F < F$. However, according to the following proposition, an admissible density 
functional Fo(p) can differ from F(p) only when d~Fo(p ) = 0: 

Proposition 5. If Fg : X —i > R U {+oo} is an admissible functional and F : X —>• R U {+oo} the Lieb 
functional, then 


- v £ d F 0 {p) -*=>- F 0 (p) = F(p) A —v £ d F(p). (28) 

Proof. Let us assume that —v £ d~F 0 (p). From Eq. (l27l) . we then have E(v) = F 0 (p) + (v \ p). Substituting 
this result into Eq. flUD, we then find that F(p) = F 0 (p) and — v £ d F(p). Conversely, if F(p) = F 0 (p) and 
— v £ d~F(p) hold, then we have E(v) = F 0 (p) + (v | p) by Eq. (|M|l . from which — v £ d~F 0 (p) follows by 
Eq. (HID. ' □ 


C. Constrained-search density functionals 

So far, we have not established a connection between the minimizing densities in the Hohenberg-Kohn 
variation principle in Eq. (O and the minimizing wave functions in the Rayleigh Ritz variation principle 
in Eq. ©■ For instance, we have not shown that each p £ d + E(v) is the ground-state density associated with 
some wave function if. Introducing a generic (abstract) constrained-search density functional, we consider in 
this section the relation between minimizers in the Hohenberg-Kohn and Rayleigh Ritz variation principles. 

Let S' be a set with elements s called states and let d : S —>• A' be a map from S to the Banach space of 
densities X. The image 


p s = d(s) (29) 

is called the density of s. However, no physical meaning is attached to s £ S or to p £ X at this point; they 
are mathematical objects- for example, p may be a (current) density or a reduced density matrix, while s 
may be a pure state or a density operator. As before, X' is the dual space of X , meaning that v £ X' if 
and only if v : X —» R is a linear and continuous operator. Thus, the pairing (v | p) is separately linear and 

continuous in v and p. We call v a potential but again no physical meaning is attached at present. 

Assume next that, for every v £ X ', we are given a map £ v : S —> R of the form 

£ v (s) = £ 0 (s) + (v | p s ), (30) 

where £q : S —> R is bounded below, meaning that there exists anMgl such that 

£ 0 (s) >M, Vs £ S. (31) 

We call £ v (s) the expectation value of the total energy when the system is in the state s and influenced by 
the potential v. This setting covers all standard formulations of DFT and CDFT, including both pure-state 
and density-matrix formulations- in particular, we do not assume linearity of £o(s) and d{s) in s. 

Next, define the ground-state energy E : X' —> R U {—oo} by the Rayleigh- Ritz variation principle: 

E(v) := inf £ v (s) = inf (£ 0 (s) + (v | p s )) ■ (32) 

s£S s£S 

Observing that v couples to p s only during minimization, it makes sense to define the constrained-search 
functional Fqs :X->KU {+oo} as the map 

-Fcs(p) = inf {£o(s) I s p} , (33) 


6 


which takes the value +oo whenever there is no s £ S with p s = p. Combining Eqs. (1321) and (1331) . we obtain 
the Hohenberg-Kohn variation principle 

E(v) = inf (Fqs (p) + (v | p)) ■ (34) 


This result proves: 

Proposition 6. The constrained-search functional Fes in Eq. (1331) is an admissible density functional for 
the energy E in Eq. m- 


D. Characterization of ground states 

Suppose that v £ X' is such that there exists s £ S for which the infimum in Eq. (1321) is a minimum: 

E(v) = inf £ v (s') = £ v (s), (35) 

s' 

meaning that there exists a ground state s £ argmin s , £ v (s r ) for v. Let now v' £ X' be arbitrary, and 
compute 

E(v') = inf £ v '(s') < £ v '(s) = £ 0 (s) + (v'\p s ) = £ v {s) + (t/ - v\p s ) = E(v) + {v' - v\p s ), (36) 

s'Go 

where we have used Eqs. (13U1) and (l35l) . We have thus proved the following proposition: 

Proposition 7. If s £ S is a ground state for v in Eq. then p s is a supergradient of E in Eq. at 

v: 

s £ argmin £ v (s') => p s £ d + E{v). (37) 

s' 

The following proposition establishes an important relationship between ground states and subgradients 
of Fes: 

Proposition 8. Let Fes : X —> R U {+oo} be the constrained-search functional (1331) . If s is the ground 


state for some potential v £ X' with density p s , then —v £ d Fes(ps)- 

Proof. Let s' £ S with p s > = p s ■ According to Eq. (13^1) . we then have 

E(v) = £ 0 (s) + (u|p s ) < £ 0 (s') + (v\p s ). (38) 

Subtracting (v\p s ) from both sides, taking the infimum on the right-hand side, and using the definition 
in Eq. (l33l) . we obtain 

£o(s ) < inf £ 0 {s') = Fcs(p s ) < ^o(s)- (39) 

S't~>p s 

Thus, if there is a ground state s £ S for v £ X\ then 

Fcs(Ps) + (v\p s ) = £ 0 (s) + (v\p s ) = £ v {s) = E(v) < F C s(p) + (v\p), VpSl. (40) 

Rearranging, we obtain 

Fcs(p) > Fcs(Ps) - {v\p - p s ), \/p£X, (41) 

demonstrating that — v £ d~Fcs{ps) and completing the proof. □ 


Conversely, does the subgradient relation —v £ d Fcs(p) imply that there exists a ground state s ^ pi 
To answer this question affirmatively, we must assume that Fes is expectation valued , defined as follows: 
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Definition 1 (Expectation-valued constrained-search functional). A constrained-search functional Fqs is 
called expectation valued if, for every p with Fcs(p) < +oo, there exists an s p £ S such that 

Fes(p) = inf £ 0 (s) = £o(s P ), s p p. (42) 

Si —V p 

In other words, if Fes is expectation valued, then the infimum in Eq. (j33[) is a minimum if Fcs{p) < +oo, 
implying that Fcs(p) is the expectation value £q (s) of some state s M- p. If Fqs is expectation valued, it 
follows immediately that 

- v£d~F C s(p ) => E(v) = F C s{p) + (v\p) => E(v) = £ 0 (s p ) + (v\p) = £ v (s p ) => s p £ argmm£ v (s') 

( 43) 

and therefore that there exists a ground state s p of v. Summarizing, we have proved the following: 

Proposition 9. Suppose that Fes • X —> MU {+oo} is an expectation-valued constrained-search functional. 
Then, for each v £ X', we have 

s p £ argmin^(s') •<=>• —v £ d~F qs{p) => p £ d + E(v). (44) 

s' 

Remark: Proposition [9] for an expectation-valued constrained-search functional should be compared with 
Proposition []] for a general admissible density functional. Whereas, for a general admissible density func¬ 
tional Fes, the condition E(v) = Fcs(p) + M/°) does not imply the existence of a ground state s 4 /), 
this implication does hold for for an expectation-valued constrained-search functional. Compare also with 
Proposition [3j where the admissible density functional is assumed to be convex and lower semi-continuous 
(i.e., the Lieb functional). 

Proposition 10 . Suppose that an expectation-valued constrained-search functional Fes ■ X —> RU{+oo} is 
convex and lower semi-continuous so that Fes = E. Then, for each v £ X’, we have 

s p £ argmin^(s') -£=>■ —v £ d~Fcs(p) -£=>■ P £ d + E[v). (45) 

s' 

Proof. Combine Propositions [3] and [HI □ 


E. Differentiability of E(v) 

We consider here conditions for differentiability of E : X ' —> R U {—oo}. From Eq. (1331) , it follows that E 
is a concave and upper semi-continuous function. Upper semi-continuity is, however, a rather weak property. 
What is needed to make E continuous ? It is a fact that a concave (convex) map over a Banach space is 
continuous on the interior of its domain Q. Hence, if we can show that E{v) is finite for every v £ X ', 
then the domain dom(E) = X' and continuity follows. The trick to show finiteness of E(v) in the context of 
standard DFT is to observe that s V„(s) = {v\p s ) is relatively bounded with £o~bound smaller than one, 
meaning that, for each v £ X', there exist e £ (0,1) and C e > 0 such that 

|V„(s)| < e£ 0 (s) + C e , Vx£S. (46) 

Assuming that V v is relatively bounded with e < 1 and with no assumption on C e , we obtain: 

£ v (s) = £ 0 (s) + V v (s) > £ 0 (s) - e£ 0 (s) - C e = (1 - e)£ 0 (s) - C e . (47) 

Taking the infimum over s £ S and using the fact that £q(s) is by definition below bounded, we find 

E{y) > (1 - e)M -C £ > -oo. (48) 

Proposition 11. Let £q : S —I>• M, E : X' —> R U {—oo}, and V v : S —> R be as described above. If V v is 
relatively bounded with respect to £q with bound e < 1 for each v £ X’, then E is continuous on X’. 


Proof. E is everywhere finite, so dom(F) = X'. Therefore E is continuous at any point v £ X'. □ 

Having determined sufficient conditions for continuity of E, what about differentiability? By definition, 
E is superdifferentiable at v £ X' if the superdifferential d + E(v) is non-empty. Clearly, if E is differentiable 
at v, the superdifferential is a singleton, d + E(v) = (VF(ii)}. The converse is in general almost true Q: 

Theorem 1. Let X be a Banach space. A convex (concave) map f : X —» lU {— oo,+oo} is Gateaux 
differentiable at x £ X if and only if it is continuous at x with a unique subgradient (supergradient). 

Hence a unique sub- or supergradient does not by itself guarantee continuity and therefore not (Gateaux) 
differentiability: A subdifferential may be a singleton even if the function is non-differentiable. This subtle 
point illustrates the limitations of finite-dimensional intuition and has been the source of misunderstanding 
in the DFT literature. For a discussion, see Ref. 0 and references therein. 

We can finally establish a useful statement on the differentiability of E : X' —>• R U {— 00} : 

Proposition 12. Let E : X' —> R be everywhere finite and given by the Rayleigh-Ritz principle in Eq. H32\) 
and assume that Fes Eq. 133\) is expectation valued and equal to the Lieb functional (i.e., convex and 
lower semi-continuous). Then E is differentiable at v if and only if all ground states s £ argmin s , £S £ v (s r ) 
supported by v have the same ground-state density, s e- > p £ X. In particular, E is differentiable at v if there 
exists a unique ground state s for v. 

Proof. Being continuous, E is differentiable at v if and only if d + E{v) is a singleton. But by Proposition HOl 
d + E{y) is a singleton if and only if all ground states of v have the same density, using the assumption that 
Fes is the Lieb functional. □ 

Remark: The message here is that given the Rayleigh- Ritz variation principle, it is not sufficient to know 
that a potential v £ X' has a unique ground-state density p = p s to prove differentiability of E at v. To 
identify the unique supergradient p £ d + E{y) with a ground-state density, Fes must be identical with the 
Lieb functional (convex and lower semi-continuous) and expectation valued. If the Lieb functional happens 
to be different from every possible constrained-search functional, differentiability may fail, even if v has a 
unique ground-state density p s . A counterexample is given in Sec. IIIIl where we consider the Hilbert space 
Xu = F 2 (R 3 ) rather than the usual Banach space Xl = F^R 3 ) fl F 3 (R 3 ). Unlike E : Xl —t R, the function 
E : Xu —» R is not differentiable. 

Starting from a Hohenberg-Kohn variation principle with an arbitrary admissible density functional, it 
seems hard to prove differentiability from the sole assumption that v has a unique ground-state density p s . 
Identifying the proper constrained-search functional, based on the Rayleigh Ritz variation principle, seems 
the only way out of the problem. 


III. APPLICATION TO DFT FOR MOLECULAR SYSTEMS 
A. Ground states and ground-state densities 

For A-electron systems, each normalized wave function has a density p £ F 1 (R 3 ) with p > 0 almost 
everywhere and fp( r) dr = N. The set of states S can be taken to be either the L 2 normalized states if £ 
Fg(R 3Ar ) (pure states) or the set of density matrices T = \k\f>k){ifk\ constructed as convex combinations 
from an orthonormal set {ifk} C Hg(R 3Ar ) with A*, > 0 and ]T) fe Afc = 1 (mixed states). In the pure-state 
and mixed-state cases, respectively, constrained search gives the functionals Fll defined in Eq. © and F D m 
defined in Eq. 

All (p) = mf £ 0 (if), £o(V0 = {i>\T + W\if), V v {ip) = {if\V\tf} = (p^\v), (49) 

Wi-Ap 

Fdm(p) = inf £b(r), £o(T) = Tr[(f + W)T], V„(T) = Tr(rU) = (prM, (50) 

r^-p 

where V is the multiplication operator associated with v. It is obvious that Fdm < Fll- It is less obvious, 
but true, that there are p £ L 1 for which Fdm(p) < Fll(p). 
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In a classic paper pjl], Kato demonstrated that, for each v £ F 2 (R 3 )+F°°(R 3 ), the X-electron Hamiltonian 
H = T + W + V is self-adjoint on L 2 (R 3Ar ) with domain U|(R 3JV ) and bounded below. More generally, 
Simon demonstrated that, for each v £ F 3 / 2 (R 3 ) + £°°(R 3 ) (which includes all Coulomb potentials), V v is 
relatively £o-bounded with an arbitrarily small bound [141 ] : 

Theorem 2. Let v £ F 3 / 2 (R 3 ) + F°°(R 3 ) and let V be the corresponding multiplication operator. Then, for 
any e > 0, there exists a constant C e > 0, such that for all if £ Ug(R 3JV ) with ||V’||2 = 1, 

mvm <e(if\T + W\if)+C e . (51) 

Proof. See Ref. [3]. □ 

In the pure-state case, therefore, each v £ F 3//2 (R 3 ) + F°°(R 3 ) gives a V v relatively bounded by So, 
with bound e arbitrarily small. Lieb proved that Fll and Fdm are both expectation valued Q. In the same 
publication, a proof (due to Simon) that Fdm : T 1 (R 3 ) —>• RU{+oo} is convex and lower semi-continuous was 
given. Using the Sobolev embedding theorem, it can be shown that for each if £ Hg(M. 3N ), the corresponding 
density p £h i/’ belongs to F 3 (R 3 ) [4j. Thus, it is appropriate to consider the Banach space 

X L := F^R 3 ) rU 3 (l 3 ) (52) 

with norm || ■ || = || • || L i + |j ■ || L 3 , whose dual space is 

Xf = F 3/2 (R 3 ) + F°°(R 3 ) (53) 

with the topology induced by || • || [l5|. Since convergence in Xl implies convergence in L 1 , Fdm ■ Xl —> 
RU {+cx)} is lower semi-continuous as well. Applying Proposition [lO] for Fdm (which is expectation valued 
and lower semi-continuous convex) and Proposition [9] for Fll (which is expectation valued but not lower 
semi-continuous convex), we obtain 

Corollary 1. Let v £ XL = L 3 / 2 (R 3 ) + L°°(R 3 ). Then, 


T £ argroin (fo(r') + (pr'l v)') 
if £ argmin (£o{4>') + v)j 


-v £ d F D M(pr) pr £ d + E(v), 
-v £ 9 _ Fll(/9j/j) => pp £ d + E{v). 


(54) 

(55) 


Note that p £ d + E(v) does not imply the existence of a ground-state wave function for v such that if p, 
only the existence of a mixed ground state T <—> p. For example, if v £ Xf has a two-fold degeneracy with 
pure ground states ifi and ip 2 , then Xpi + (1 — X)p 2 £ d + E(v) but there may be no pure ground state with 
this density. On the other hand, X\ifi)(ifi \ + (1 — X)\if2)(if2\ is a mixed ground state with this density. 


B. Topology dependence of the Lieb functional 

Recall that the Lieb functional F :X ^RU {+°o} depends explicitly on the Banach space X, including 
its topology. We now demonstrate this using an explicit example. Motivated by the set inclusion 

X L = F 1 (R 3 ) n F 3 (R 3 ) C X h := L 2 ( R 3 ) (56) 

and the fact that Xh is Hilbert space with a simpler structure than the non-reflexive space Xl, it is tempting 
to consider Fdm as a function on X = Xp. However, the embedding Xl C Xp is not continuous: Conver¬ 
gence in Xl does not imply convergence in Xp. Even if Fdm is lower semi-continuous with respect to Xl, 
it may fail to be so in Xp. To see this, we note that, since Xp is a Hilbert space, X(j = Xp, which does not 
admit constant potentials: If v £ Xp, then v + c ^ Xp for each constant 0. This observation allows us 
to prove the following result: 
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Proposition 13. For each v £ Xh, we have E{y) <0. If v > 0 almost everywhere, then E{v) = 0. 
Proof. See Appendix A. 


□ 


The proof is based on the following idea. Write v = — V- where > 0 almost everywhere. The 

negative part V- can only lower the energy, whereas v+ £ Xh implies that v decays at infinity (in an average 
sense), thereby allowing the electrons to lower their energy to zero by “escaping to infinity”. Note how the 
fact that nonzero c ^ Xh affects this argument. 

Consider now p = 0, for which Edm(p) = +oo since no T i-> 0. Evaluating the conjugate (Lieb) functional 
with respect to X h at p = 0, we obtain 

F(0) = sup (E(v) — (u|0)) = sup E{y) = 0. (57) 

v£X h v£X h 

Note that a different result is obtained using Xl, which admits constant potentials v(r) = c £ R: 

F( 0) = sup (E(v) — (u|0)) = sup E(v) = +oo = Edm(0). (58) 

v£X l v£X l 


As an immediate consequence, we arrive at the following result on Xh: 

Proposition 14. For any v £ Xh such that v > 0 almost everywhere, it holds that, 

E{v) = F{ 0) + (u|0) A —v £ d~F{ 0) A 0 £ d+E{v), (59) 

where F is conjugate to E. There exists no ground state i/> £ Hg(M. 3N ) for any v > 0. 

Proof. Eq. (l59l) follows from Proposition[3] It remains to show that there exists no ground state if v > 0 
almost everywhere. Since E(v) = 0, it is sufficient to show that, for each if £ iEg(R 3JV ), it holds that 
£ v {ip) > 0. We have £ v (ip) > (i/j\T\ip). But if if has zero kinetic energy, then Vif = 0, from which it follows 
that if = 0 almost everywhere, contradicting HV'lh = 1- Thus £ v (ip) > (i/>|T|i/>} >0. □ 

Our discussion illustrates how the choice of Banach space X of densities affects the properties of the 
conjugate universal functional F : X — > R U {+oo}. Clearly, the Lieb functional Edm for the space Xl is 
very different from the Lieb functional for Xh, since, in the latter case, for any non-negative potential v, the 
unphysical density p = 0 is a minimizer of the Hohenberg-Kohn variation principle. But v does not even 
have a ground state and certainly not a ground state with density p = 0. 

Connecting with the discussion following Proposition ! 121 we can see how topology influences differentia¬ 
bility of E. The map E : Xh -a R is pointwise identical to E : Xf —» R: it is by definition the ground-state 
energy of the system in the external potential v, defined via the Rayleigh Ritz variation principle. However, 
the topology on Xh is different and Fdm is no longer equal to the Lieb functional. We then cannot conclude 
from PropositionfTSl that E : Xh -a R is differentiable at v if v £ Xh has a unique ground state. 


IV. APPLICATION TO CDFT FOR MOLECULAR SYSTEMS IN MAGNETIC FIELDS 

Consider an X-electron system subject to an external magnetic vector potential A : R 3 —> R 3 , with 
associated magnetic field B(r) = V x A(r) (in the distributional sense). Considering A as a variable on 
the same footing as the scalar potential v, we arrive at CDFT, where the paramagnetic current density 
j p £ L 1 (R 3 ) appears as a variable together with p [lj, [l6}. The mathematical foundation of CDFT is 
not as well developed as that of DFT. Indeed, part of the motivation for the present work stems from a 
study of CDFT [lj. On the other hand, Laestadius [§} has taken important steps—proving, for example, 
Proposition [T7] below. 

Following Ref. j§], we assume for simplicity that the components of A are L°°(R 3 ) functions. The single¬ 
electron momentum operator —iV is then replaced by —iV + A, transforming the X-electron kinetic-energy 
operator T into the corresponding Ta- Since A £ L°°(R 3 ), the kinetic-energy expectation value is still well 
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defined for each if £ Hg(M. 3N ). In terms of the density p £ L 1 (K 3 ) and the paramagnetic current density 
j p £ L 1 (R 3 ), we obtain 

(0|Ta|0) = (0|T|0) + \ j I A(r)| 2 p(r) dr + j A(r) • j p (r) dr. (60) 

The ground-state energy is thus given by 

E {v, A ) = , inf , (0 I^a + W + V\if>) = inf (F VR (p, j p ) + (v + \A 2 \ p) + (A | j p )) , ( 61 ) 

where we have introduced the Vignale-Rasolt (VR) functional (§] 

^Vr(pJ p ) := inf |(0|T + W\if>) j 0 £ H^(R 3N ), [|0|| 2 = 1, 0 ^ (p,j P )} • (62) 

Similarly, we may define a density-matrix constrained-search functional 

Fum(p, jp) = inf {Tr((T + W)T) \ T h- (p, j p )} , (63) 

where pr = YhkPkPi’k anc i Jpr = Sfe Pkjp,j, k • Whereas Fdm is convex by construction, the presence of the 
nonlinear A-dependent term makes E{v,A) nonconcave. Following Ref. fl6| . it is therefore natural instead 
to work with E(u , A) = E(u — \A 2 , A), defined in Proposition [TTjl below. We note that 

(p,jp)£A L xL 1 (R 3 ), (64) 

a Banach space with norm ||(p, j p )|| = ||p|| + ||j p H l 1 and dual x L°°(R 3 ), the latter containing all (v , A). 
First, we prove finiteness of the ground-state energy: 

Proposition 15. For every (v,A) £ X R x L°°, the ground-state energy is finite, E{v, A) > —oo. 

Proof. For an A-electron system influenced by a magnetic field, the diamagnetic inequality [l8j | gives 

(01 T\<f) < (0|Ta| 0), V0 e H'(R 3N ), (65) 

where 0 = |0| be the pointwise absolute value. Next, using the identity p^ = p$ and the relative boundedness 
of v £ X R with respect to T, we find that V is relatively bounded with respect to Ta: 

1(01^10)1 = 1(01^10)1 < e(0|f |0) + < c(0|Ta|0) + C e . (66) 

Similarly, W is relatively bounded with respect to Ta- It follows that (0|.Hj0) is below bounded, and thus 

that E(v, A) > —oo. □ 

Next, we note that, for each A £ L°°(R 3 ), we have |A| 2 £ L°°{\ R 3 ) C X[ . Therefore, for each pair 
(v, A) £ X R , we have (v ± I|A| 2 , A ) G X L x L°°, making the following proposition easy to prove: 

Proposition 16. If F 0 : Xf x L 1 is either Fy R or Fq m, then E : X £ x L°°(R 3 ) —>• R defined by 

F{u, A) = inf (Fo(pJp) + (« I p) + (A|j p )) (67) 

(p,Jj>) 

is concave, finite, and therefore continuous. The ground-state energy E and E are related by 

E(v,A)=E{v + ±\A\ 2 ,A), E(u, A) = E(u — I|A| 2 , A). (68) 

Proof. We leave the details to the reader. □ 
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The map (v, A) (v + i | A| 2 , A) defined on X £ x L°°(R 3 ) is smooth and invertible, with smooth inverse 
{u, A) i-> (u — i | A| 2 , A). Thus, the properties of E are reflected in properties E and vice versa. If F 0 is an 
admissible functional for E : X £ x L°° — > R, then 

E{v, A) = Fo(p, j p ) + (v + i|A| 2 |p) + (A|j p ) •<=>■ — (v + ||A| 2 , A) e d i'b(pjp) ^ 

(pJp) S d + E(v + ^|A|", A) 

Laestadius proved the following result js|: 

Proposition 17. FVr : T 1 (R 3 ) x LfyR 3 ) ->KU {+ 00 } is expectation valued. 

Proposition [9] therefore gives 

G argmin (ip\H\ij}) E(v, A) = F V r(pJ p ) + (v + \\M 2 \p) + ( A |j P ). (70) 

i>' 

Thus, using the Vignale-Rasolt functional and the space Xl x L 1 as density space, there are no minimizers 
of the pure-state Hohenberg-Kohn variation principle for CDFT that do not correspond to ground-state 
wave functions. 

On the other hand, it is not known whether the density-matrix functional Tdm in (162H is expectation 
valued, but it seems likely. If Tdm is not expectation valued, there may be (p, j p ) such that 

E{v, A) = Fdm(Pj j P ) + (v + ^|A| 2 |p) + (A|j p ) (71) 

but such that there is no T i-> (p, j p ), i.e., (p, j p ) must be considered an unphysical minimizer of the density- 
matrix Hohenberg-Kohn variation principle for CDFT. 

Neither is it known whether Fdm is lower semi-continuous in the L 1 x L 1 topology; the proof for the 
ordinary DFT case in Ref. 0 is not easy to generalize to the present situation. Thus, there could be 
(p,j p ) € d + E(v + ^|A| 2 ,A) such that — (v + ^|A| 2 ,A) ^ <9 - Fdm(p, j P ), so that (p,j p ) does not satisfy 
Eq. (1711) . 

Finally, we may consider the Lieb functional F : Xl x L 1 —>• RU {+ 00 }, defined as 

F(P Jp) := sup (e(u, A) - (u\p) - (A|j p )) = sup (E(v, A) - (v + i|A| 2 |p) - (A|j p )) . (72) 

This functional is convex and lower semi-continuous by construction. However, unlike in standard DFT, it 
is unknown whether F(p, j p ) = F D m(p, jp)- 

From the perspective of applying Proposition [T01 thereby establishing a result like Corollary [[] for CDFT, 
one must establish both that Fdm is expectation valued and lower semi-continuous. Thus, we do not know at 
the present time, whether E(y, A) is Gateaux-differentiable when the ground-state density (p,j P ) is unique. 


V. CONCLUSION 

We have analyzed the relationship between ground-state densities obtained from the Rayleigh Ritz varia¬ 
tion principle (by minimizing over pure-state wave functions or mixed-state density matrices) and from the 
Hohenberg-Kohn variation principles (by minimizing over densities using an admissible density functional 
Fo(p)). For standard DFT for molecular systems, we established, for each potential v £ L 3 / 2 (R 3 ) +L°°(R 3 ), 
a one-to-one correspondence between the mixed ground-state densities of the Rayleigh-Ritz variation princi¬ 
ple and the mixed ground-state densities of the Hohenberg-Kohn variation principle with the Lieb density- 
matrix constrained-search functional Fdm : L 1 (R 3 ) D L 3 (R 3 ) ->1U {+ 00 }. A similar one-to-one correspon¬ 
dence is established between the pure ground-state densities of the Rayleigh Ritz variation principle and 
the pure ground-state densities of the Hohenberg-Kohn variation principle with the Levy-Lieb functional 
Fll : L 3 (R 3 ) D L 3 (R 3 ) —>• R. In other words, all physical ground-state densities (pure or mixed) are recov¬ 
ered with these functionals and there are no false densities (i.e., minimizing densities not associated with a 
ground-state wave function). 


13 


We also noted how the topology of the underlying Banach space X impinges on the results— in particular, 
we noted that the Lieb functional F : X —> RU{+oo} depends explicitly on X. As an illustration, F ^ Fdm 
on X = L 2 (R 3 ) but F = Fdm on L 1 (M 3 ) n L 3 (K 3 ). Finally, CDFT was discussed and some open problems 
were pointed out—for example, it is unknown whether FbM(p,jp) is lower semi-continuous in any useful 
topology such as L 1 (R 3 ) x L 1 ' 
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Appendix A: Proof of Proposition 1131 


Proof of Proposition Wft Writing v = — V -, we obtain 

E(v) = inf (Fdm(p) + ( v+ \, p) - (u_ |,p)) < inf (F DU (p) + (v+ \ p) = E(v +)). (Al) 

p p 

It is therefore sufficient to show that E(v+) < 0. In fact, we show that E(v+) = 0. 

Let v > 0 almost everywhere. Let A > 0 be arbitrary and let H^a with k £ N be disjoint cubes of side 
length A such that = R 3 - Each can be obtained by translation of H^a- Since v £ L 2 (R 3 ), 

/ u(r ) 2 dr = / v(v) 2 dr < +oo, 

J& k,\ 

implying that f Qk x v(r) 2 dr —> 0 as k — > oo. Let ^ £ Cf°(M. 3N ) be arbitrary but with support in Ufj so that 
PtI, £ C(?°(R 3 ) has support contained in By translating if properly (denoting the result by ijj k ), the 

support of ifk is inside it kl and 

Edm(pvJ = Fdm(PiIi) < (i/>\T + W\il>) = ( T) + (W), (A2) 

independent of k. We obtain 


E(v) < lim ^(T) + (I V) + v(r) p k { r) drj = (T) + (W), 


where we have used the fact that 


/ v(r) Pfe(r) dr < ( / u(r ) 2 dr I 


1/2 


-A 0 


(A3) 


(A4) 


as k —> oo. 

We now increase the size of the boxes Q k} \ by varying A > 0. By dilating if in the manner 

^(ri,---) ^ A 3Ar/2 i/>(Ari,---), (A5) 

the support is still inside and the density is scaled as p^(r) -A A _ 3 p^(A _ 1 r), conserving the number of 
particles. We obtain the scaling 

(T) + ( W) -A \~ 2 {T) + A _ 1 (VF). (A 6 ) 


By repeating the above argument for A = 1 and letting A —> oo, we obtain E{v) < 0. On the other hand, 
E{v) > 0 since the Hamiltonian H{v ) is positive with v > 0, yielding E(v) =0. □ 
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